const int Nch = 5; 
double BunchLength = -1.;
int Target = -1;
double tmin_beam_ped = 0.; double tmax_beam_ped = 0.;
bool SelCh[Nch] = {false, false, false, false, false};
double tmin_beam_sgn = 0.; double tmax_beam_sgn = 0.;
double tmin_bkg, tmax_bkg;
double Ti0,Tf0,Ti1,Tf1,Ti2,Tf2;
bool Alignment = false;
int ich_ref;
double T_al_1, T_al_2;
int ReferenceSignal; 
double tmin_plot,tmax_plot;
double Ne_norm;
int  FreqCutType = 0; 
double FreqCutMin = 2;   
double FreqCutMax = 6;  


double Ti_all=-306 * ns, Tf_all=-296.5 * ns;
double Fmin_av = 0.001, Fmax_av = 2.;
double Ne_min = 0., Ne_max = 0.;
bool DoCalibration;
TGraph *tgCalib;
TGraph *tgAntenna;

double FrascatidEdX[7] = {1.25303e+09, 
			  1.56448e+09, 
			  1.87671e+09, 
			  1.29787e+09, 
			  7.90366e+08, 
			  4.0957e+08, 
			  2.22723e+08}; //[MeV], Ne=1e9


double SLACdEdX[7] = {4.32413e+08, 
		      8.68644e+08, 
		      3.66211e+09, 
		      7.8667e+09, 
		      1.19264e+10, 
		      1.44541e+10, 
		      1.3692e+10}; //[MeV], Ne=1e9




double FrascatidEdXFoV[7] = {4.69284e+08, 
			     5.6268e+08, 
			     6.09154e+08, 
			     3.9658e+08, 
			     2.31763e+08, 
			     1.20065e+08, 
			     6.4932e+07}; //[MeV], Ne=1e9

double SLACdEdXFoV[7] = {2.88262e+08, 
			 5.67987e+08, 
			 2.30111e+09, 
			 4.80265e+09, 
			 7.09747e+09, 
			 8.38416e+09, 
			 7.77865e+09}; //[MeV], Ne=1e9



double FrascatiStoppingPowerCenter[7] = {108.225, 
				135.27, 
				157.501, 
				106.97, 
				63.479, 
				33.1579, 
				17.4653}; //[TeV cm2 / g], Ne=1e9

double SLACStoppingPowerCenter[7] = {90.725, 
				     391.829, 
				     1625.8, 
				     3492.12, 
				     5251.11, 
				     6395.3, 
				     6062.93}; //[TeV cm2 / g], Ne=1e9

/*
double FrascatidEdXCenter[7] = {220.337, 
				108.582, 
				32.2893, 
				12.6924, 
				7.32219, 
				4.5989, 
				0.933167}; //[MeV], Ne=1e9

double SLACdEdXCenter[7] = {241.014, 
			    439.468, 
			    1360.97, 
			    2298.91, 
			    2818.21, 
			    2738.93, 
			    2161.87}; //[MeV], Ne=1e9
*/

//1 dm3
double FrascatidEdXCenter[7] = {49.0527, 
				24.4837, 
				6.57785, 
				2.96432, 
				1.57843, 
				0.678571, 
				0.304067}; //[MeV], Ne=1e9

double SLACdEdXCenter[7] = {80.304, 
			    146.221, 
			    451.853, 
			    761.92, 
			    931.852, 
			    903.372, 
			    711.787}; //[MeV], Ne=1e9





 


